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We calculate the nucleon strange quark content directly from disconnected three-point functions. 
Numerical simulations are carried out in two-flavor QCD using the overlap quark action with 
up and down quark masses down to a fifth of the physical strange quark mass. To improve the 
statistical accuracy, we calculate the nucleon two-point functions with the low-mode averaging 
technique, whereas the all-to-all quark propagator is used for the disconnected quark loop. We 
obtain the y parameter, which is the ratio of the strange and light quark contents, y = 0.024(45) 
at the physical point. This is in a good agreement with our earlier calculation from the nucleon 
spectrum through the Feynman-Hellmann theorem. 
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1. Introduction 



The nucleon strange quark content (N\ss\N) is an important parameter to determine the cross 
section of the scattering of dark matter candidates from the nucleon [jl|, ||]. It can not be measured 
directly by experiments, and only lattice QCD can provide a model-independent and nonpertur- 
bative determination. A precise lattice calculation is, however, very challenging, because only 
disconnected diagrams contribute to the strange quark content and they are computationally very 
expensive to calculate with the conventional method. In addition, the scalar operator ss has a vac- 
uum expectation value (VEV), which diverges towards the continuum limit. We need to subtract 
the VEV contribution and this induces a substantial uncertainty in the strange quark content. 

In our previous study we avoid the above mentioned difficulties by calculating the strange 
quark content from the quark mass dependence of the nucleon mass through the Feynman- 
Hellmann theorem 

(N\ss\N) = (1.1) 
om s 

We refer to this method as the spectrum method in the following. This method is, however, not 
applicable to other interesting matrix elements, such as the strange quark spin fraction of the nu- 
cleon. In this article, therefore, we attempt a direct determination of strange quark content from 
nucleon matrix element including a disconnected diagram. To this end, we employ the overlap 
quark action, which has exact chiral symmetry, and improved measurement methods, such as the 
low-mode averaging (LMA) technique [Q, |5|] and the use of the all-to-all quark propagator 



2. Simulation details 

Gauge ensembles of two-flavor QCD are generated on a L 3 x T = 16 3 x 32 lattice using the 
Iwasaki gauge action and the overlap quark action. We set the gauge coupling j8 = 2.30 at which the 
lattice spacing determined from the Sommer scale r$ = 0.49 fm is a = 0.1 18(2) fm. Our simulation 
is accelerated by introducing a topology fixing term into our lattice action [0], and we simulate 
only the trivial topological sector Q = in this study. We take four values of bare up and down 
quark masses m uc i = 0.015,0.025,0.035 and 0.050, which cover a range of the pion mass m K = 
290 — 520 MeV Statistics are 100 independent configurations at each quark mass. We refer readers 
to [^] for further details on our configuration generation. In our measurement, we take two values 
of the valence strange quark mass m s va i = 0.070 and 0.100, which are close to the physical mass 
ms,phys = 0.077 determined from our analysis of the meson spectrum 

The strange quark content can be extracted from nucleon two- and three-point functions 

C Y 2pt {tAt) = Tv[T(N{t + At)N{t))\, (2.1) 
C$ pt (t,At,At s ) = Tr[r(Ar(? + Af)5(f + A/ s )iV(/)}]-(5(/ + A/ s ))Tr[r(7V(f+Af)^(?)}], (2.2) 

where S = ss is the strange scalar operator, t represents the temporal coordinate of the nucleon 
source operator, and At (At s ) is the temporal separation between the nucleon source and sink (quark 
loop). We calculate t and C\ pt with two choices of the projector r = F± = (1 ± j4.)/2 corre- 
sponding to the forward and backward propagation of the nucleon. We then take the average over 
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the two choices of T with appropriately chosen temporal separations At and At s . The averaged cor- 
relators, which are denoted by C2 P t and C^ pt in the following, show reduced statistical fluctuation. 

For further improvement of the statistical accuracy, we employ the low-mode averaging(LMA) 
technique [Q, ||] to calculate C2 P t and C^, pt . In this method, the quark propagator is expanded in 
terms of the eigenmodes of the Dirac operator D. We calculate the contribution of 100 low-modes 
exactly as 

100 i 

(°)r: = I^v«v«t, Dv«=A«v«. (2.3) 

The remaining contribution from the higher modes is taken from by that of the conventional point- 
to-all propagator. With this decomposition of the quark propagator, C2 pt is divided into eight con- 
tributions 

<-2pf — + c 2p/ ~r^2 P t + c 2pf ~r c 2pf + c 2p; + c 2pf + ( -2p/ • \ L -^) 

It is expected that C^f dominates C2 P t at large temporal separation At. The statistical accuracy 
of C2 P t can be remarkably improved by averaging C^f over the location of the nucleon source 
operator. We also improve the statistical accuracy of other contributions (Cl^f, C™ 77 ) by using 
point-to-all propagators averaged over 4 or 8 different source locations. The nucleon piece of the 
disconnected correlator Cj, pt is calculated in the same way. 

Since the nucleon correlators C2 P t and Ci, pt damp rapidly as At increases, it is essential to 
reduce the contamination from excited states at small At. In this study, we employ the Gaussian 
smearing 

f / CO \ N> \ 3 
q sm rM=L{( l + ^ H ) ) ?(y> f )> tfx,y = XJ(S x , y - 1 + S Xiy+ ,-) (2.5) 

for both of the source and sink operators. The parameters ft) = 20 and ,/V = 400 are chosen so that 
the effective mass of C2 P t shows a good plateau. For comparison, we repeat our measurement with 
the local sink operator. In this additional measurement, we test the local and an exponential source 
operator q smr (x,t) = £r ex P( — B\r\)q(x + r,t). The parameter B is chosen so that the distribution 



of the smeared quark is close to that of the Gaussian smearing (2.5). 

To calculate the disconnected quark loop in C^ pt , we construct the all-to-all quark propagator 



as proposed in |QJ. The low-mode contribution is the same as in (2.3) and the contribution from 
the high-modes is estimated by employing the noise method with the dilution technique We 
prepare a single Z2 noise vector v\ for each configuration and it is split into = 3x4x7/2 
vectors 7] w (d = l,...,Nd) which have non-zero elements only for single color and spinor indices 
and two consecutive time-slices. The high-mode contribution is then given by 

(»)^ = IV^V^, (2-6) 

d=\ 

where is the solution of the linear equation 

D^ = (l-^)i|M (2.7) 
and £?>i ow is the projection operator to the subspace spanned by the low-modes. 
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Figure 1: Ratio R(At , At s ) at m u( i = 0.050 and m s va i = 0.100 with Af fixed to 10. Triangles show data with 
Gaussian smeared source and sink, whereas circles (squares) are with the local sink and local (exponentially 
smeared) source. We omit the noisy high-mode contribution to the quark loop in this plot. Vertical lines 
show the location of the source and sink operators. 



LLH 



We also tested the all-to-all quark propagator to calculate the high-mode contributions C^/ 



r HHH 



in (2.4). It turned out, however, that these contributions have large statistical error due to the in- 
sufficient number of the noise samples. We therefore use C2 P t and Cj, pt calculated with the LMA in 
the following analysis. 

3. Matrix element at simulated quark masses 

We extract the unrenormalized matrix element (N\ss\N)i at from the ratio 

C 3pt (t,At,At s ) 



R(At,Al s ) 



A/,A/, 



(N\ss\N) lat . 



(3.1) 



2 P t(t,At) 

Figure [I] shows At s dependence of R(At,At s ) at our heaviest ud quark mass m ut ] =0.050 with a fixed 
At. We observe a clear plateau between the nucleon source and sink with the Gaussian smeared 
operator. On the other hand, the plateau is unclear if the local operator is used for the source and/or 
sink. It is therefore crucial for a reliable determination of (N\ss\N}iat from /?(A?,A? iS ) to reduce 
contamination from the excited states by appropriately smearing the nucleon operators. 

The situation is similar at two smaller ud quark masses m uc i = 0.035 and 0.025. As shown in 
Fig. however, we do not observe a clear signal even with the smeared source and sink at our 
smallest quark mass m U( i = 0.015. To observe a clear plateau of R(At,At s ) at such small m uc i, we 
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Figure 2: Ratio R(At,At s ) at m ud = 0.015 and m.„.„/ = 0.100. 
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Figure 3: Ratio R(At) as a function of At. Three left (right) panels show data at m s va [ = 0.070 (0.100). The 
horizontal lines show the constant fit to R(At) in terms of At to determine (N\ss\N)i at . 



may need more statistics as well as a larger lattice to suppress finite volume corrections, which are 
possibly sizable at M n L ~ 2.8 at m„d = 0.015. We leave such a calculation for future studies, and 
omit data at this m ui j in the following analysis. 

In this report, (N\ss\N)i at is determined by the following simple two-step fits. First, we carry 
out a constant fit to R(At,At s ) in terms of At s . The fit result, which we denote by R(At), is plotted 
as a function of At in Fig. [| We then extract {N\ss\N)i at by a constant fit to R(At) at 12 < At < 15. 
As seen in Fig. ||, R(At) do not show significant At dependence with this range of At. We therefore 
expect that extraction of ground state signal is well under control. 



4. Strange quark content at physical point 

As seen in Fig. ||, the fit result for (N\ss\N)i at does not have significant m s va ; dependence at 
each m u d. This leads us to interpolate (N\ss\N)i at to the physical strange quark mass m sp h ys = 0.077 
using a linear form in terms of m s , va i. Fit results as well as fitted data at m s , va i are plotted as a 
function of m ut i in Fig. |]. 



At next-to-leading order of heavy baryon chiral perturbation theory (HBChPT) [jlOJ the nu- 
cleon mass can be written as m^ = mo + C\m ui i + C2>n s + C^rn^ + C^m\ + C^m^, where Q 
(i = 1, ...,5) are functions of the low-energy constants (LECs) in HBChPT. We note that the contri- 



butions from decuplet baryons are neglected. The Feynman-Hellmann theorem (1.1) then implies 
that the m s dependence of (N\ss\N) comes from the 0{m\ ^ ) terms from K and v\ loops. By using 
the leading order relation m\ °< m u d + m s and ml oc m ut i + 2m s , we obtain 

(N\ss\N) = ^L= D +D imud + 0(m 2 u d ), (4.1) 
dm s 

where Do and Di depend on the LECs and m s . 

We extrapolate {N\ss\N)i at at m SiP h ys by the linear form ( |Q| ) with Do j treated as fitting param- 
eters. This chiral extrapolation is plotted in Fig. ^. We obtain (N\ss\N) = 0.11(21) at the physical 
point, where the error is statistical only. This is converted to the phenomenologically relevant 
parameters 

m N 
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Figure 4: Chiral extrapolation of (N\ss\N)i at at m s „hys as a function of m ui j. For a comparison, we also plot 
{N\ss\N)i ut at simulated strange quark masses m sva i =0.070 and 0.100. 



and 



2(N\ss\N) 
(N\uu + dd\N) 



0.024(45), 



(4.3) 



where we use the nucleon mass m# [11] and the ud quark content (N\uu + dd\N) obtained in our 
previous study 

As shown in Fig. ||, we observe a good agreement with our estimate from the spectrum method 
[g]. The same figure also shows that previous studies with the Wilson-type fermions []l2|, |l~3|, 14] 
led to rather large values for the strange quark content y. It is argued in |g[ |l5|] that the explicit 
chiral symmetry breaking induces a mixing between the scalar operators of sea and valence quarks 
and leads to a substantial uncertainty in the strange quark content. 
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Figure 5: Comparison of y parameter with previous estimates |^, |l2[ [TJ, |l^, |T^]. Our previous study [^J 
employs the overlap action, whereas other studies use the Wilson-type quark actions. Two results from 
the UKQCD study JTs] ] are obtained with and without subtracting the unphysical effect due to the operator 



mixing. 
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5. Conclusion 

In this article, we report on our calculation of the strange quark content directly from the 
nucleon matrix element. We determine fj t and y with an accuracy of 0(1O~ 2 ). The key points 
leading to this accuracy are the use of the improved measurement techniques, namely the LMA 
and the all-to-all quark propagator, as well as the appropriately smeared operator both for nucleon 
source and sink. It is an interesting subject in the future to extend this study to other matrix elements 
containing disconnected diagram such as the quark spin fraction of the nucleon. 

We observe a good agreement with our previous estimate from the spectrum method. Chiral 
symmetry preserved by the overlap action plays a crucial role in avoiding the unwanted operator 
mixing for the Wilson-type actions. For more precise determination, we need to extend our calcu- 
lation to Nf ■ = 2 + 1 QCD and larger volumes. Our preliminary estimate with the spectr um method 



is reported at this conference [16]. A direct determination from nucleon disconnected functions in 



Nf = 2+1 QCD is also in progress. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 09-05). This work is supported in part by the Grant-in-Aid of the Min- 
istry of Education (No. 19540286, 20105001, 20105002, 20105003, 20340047, 21674002 and 
21684013). 
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